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Abstract. Long-duration gamma-ray bursts (GRBs) provide a premier tool for studying high-redshift star-forming galaxies 
tlianks to tiieir extreme brightness and association with massive stars. Here we use GRBs to study the galaxy mass-metallicity 
(M,-Z) relation at z ~ 3 — 5, where conventional direct metallicity measurements are extremely challenging. We use the 
interstellar medium metallicities of long-duration ORB hosts derived from afterglow absorption spectroscopy (Z ~ 0.01 — 1 
Zq), in conjunction with host galaxy stellar masses determined from deep Spitzer 3.6 jlm observations of 20 ORB hosts. We 
detect about 1/4 of the hosts with Mab(^) ~ —21.5 to —22.5 mag, and place a limit of Mab (^) — 19 mag on the remaining 
hosts from a stacking analysis. Using a conservative range of mass-to-light ratios for simple stellar populations (with ages 
of 70 Myr to ~ 2 Gyr), we infer the host stellar masses and present the galaxy mass-metallicity measurements at z ~ 3 — 5 
((z) ^ 3.5). We find that the detected GRB hosts, with 2 x lO'^ Mp^, display a wide range of metallicities, but that the 
mean metallicity at this mass scale, Z ~ 0.1 Z,;;,, is lower than measurements at z < 3. Combined with stacking of the non- 
detected hosts (with < 4 X lO'' Me and Z < 0.03 Z©), we find evidence for the existence of an M-t-Z relation at z ~ 3.5 
and continued evolution of this relation to systematically lower metallicities from z ~ 2. 
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INTRODUCTION 

The simple "closed-box" model of galaxy evolution [1] predicts a correlation between the stellar mass and the gas- 
phase metallicity of a galaxy (the M*-Z relation). Such a relation is indeed seen in the nearby universe {2-A]. and 
is found to be in place out to z ^ 2 [5], progressively evolving to lower metallicities. Tracing the M*-Z relation and 
its evolution to even earlier times will provide insight into the earliest epochs of galaxy evolution, while allowing us 
to probe the relative importance of various galactic-scale phenomena proposed at z < 3. However, studying the M*-Z 
relation beyond z ^ 3 .5 is challenging because the nebular emission lines required for robust metallicity measurements 
(e.g.. Ha, Hj3, N ii A6583, iii AA4959,5007, ii AA3726,3729) shift into different near- and mid-IR wavebands, 
where existing spectrographs have reduced sensitivity compared to the optical band. This is aggravated by the faintness 
of the galaxies, making individual abundance measurements nearly impossible. 

An alternative way to determine metallicities at z > 3 (and in principle at z ^ 10 and beyond [6, 7]) is absorption 
spectroscopy of gamma-ray burst (GRB) optical/near-IR afterglows. Long-duration GRBs are known to be associated 
with the deaths of massive stars (e.g. [8]), and therefore with sites of active star formation. The large optical 
luminosities of GRB afterglows and their intrinsic featureless spectra provide a unique way to measure interstellar 
medium (ISM) metallicities for galaxies at z > 2 from rest-frame ultraviolet metal absorption lines and Ly a absorption. 
Since the afterglows are significantly brighter than the underlying host galaxies, this technique allows us to measure 
metallicities independent of the galaxy brightness. This approach has now been exploited at least to z ^ 5 using optical 
spectra (e.g., [9-11]), and with near-IR spectrographs it can be implemented to z ^ 20. 

Naturally, to explore the M*-Z relation at z ^ 3 we also require a determination of the GRB host galaxy stellar 
masses, and hence follow-up infrared observations with the Spitzer Space Telescope to probe the rest-frame optical 
luminosity. We present the first large set of Spitzer observations for GRB host galaxies at z ~ 3 — 5, and combine the 
inferred masses with measured metallicities to explore the M*-Z relation beyond z 3. A full analysis is presented in 
Laskar et al. (201 1), to be submitted to ApJ. 



GRB SAMPLE AND DATA ANALYSIS 



We obtained deep observations (about 2 hr per target) of 20 long-duration GRB host galaxies in the range z « 3 — 5.8 
using the 3.6 [im band of the Infra-Red Array Camera (IRAC; [12]) on-board the Spitzer Space Telescope. The 
effective wavelength of the IRAC 3.6 /im band probes the rest-frame spectral energy distribution (SED) redward 
of about 5500 A and therefore provides a robust measure of the stellar mass. 

We used optical afterglow images to perform relative astrometry on the Spitzer mosaics and to locate the GRB 
hosts. At the depth of our observations, Spitzer images are confusion-limited for faint sources, requiring us to model 
and subtract nearby contaminating sources in 3 cases. We detect 5 GRB hosts at 3.6 jj.m (GRBs 050319, 050814, 
060707, 060210, and 060926), with flux densities ranging from about 0.55 to 1.65 nJy, with a typical upper limit of 
about 0.25 /iJy (3ff). To assess the typical flux density of the non-detected hosts we carry out a stacking analysis of 
1 1 non-detections with accurate relative astrometry. The stacked image does not show a detection down to a 3a limit 
of < 80 nJy. 




Stellar Masses of GRB Hosts at z ~ 3 - 5 

Computing stellar masses from observed luminosities in a given wave-band requires knowledge of the mass-to- 
Ught ratio and hence the stellar population age and metallicity. When multi-band photometry is available, modeling 
of the spectral energy distribution (SED) using stellar population synthesis models can be used to determine stellar 
masses, provided that a single stellar population is assumed. When multi-band photometry is not available, the resulting 
uncertainty in the mass-to-light ratio (e.g., at ~ 1 jxm) is about an order of magnitude (e.g., [13]). 

Here, since we lack broad-band photometry, we determine a range of mass-to-light ratios for each galaxy in the 
observed 3.6 jim band using a wide range of population ages and the single stellar population models of Maraston 
[14] with a Salpeter IMF and assuming an instantaneous burst of star formation (e-folding time, T = 0). As expected, 
the 3.6 /xm mass-to-light ratio for these models increases with stellar population age beyond ^ 10 Myr An upper 
bound on the mass-to-light ratio is achieved by setting the stellar population age to the age of the universe at each host 
redshift (» 1.8 Gyr at the median redshift of our sample). We stress that this leads to a very conservative maximum 
mass for each host galaxy since studies of Lyman Break Galaxies (LBGs) and Lya emitters (LAEs) at z > 3 indicate 
typical population ages of ~ 0.1 —0.6 Gyr [13, 15-17]. The median age for long-duration GRB hosts atz ~ 1 of about 
70 Myr [18] is probably more typical. The variation in mass-to-light ratio between these age values is about an order 
of magnitude. 

The maximal masses inferred for our sample are (2.5 —5.8) x lO"' M©, while the typical (maximal) upper limits 
are < 9 X 10^ Mg . The masses inferred for a 70 Myr old population are about (0.6 — 1.4) x lO'" M©, with typical 
upper limits of < 2 x 10^ M©. The mass Umit from the stack of 11 GRB hosts is < 7 x 10^ M© for a 70 Myr old 
population, and < 3 x 10^ M© for the maximal age. We also include a literature sample of 5 long-duration GRB hosts 
at z > 3, incorporating two detections with maximal masses of 1.4 x lO^" M© (GRB 060510B; [19]) and 6.7 x lO^^ 
M© (GRB 080607; [20]). 



THE MASS-METALLICITY RELATION AT z ~ 3 - 5 



A typical optical afterglow spectrum exhibits a wide range of ISM absorption features due to rest-frame UV transitions 
of low- and high-ionization metal species, which allow a direct determination of the column density of these elements 
along the GRB line of sight through the host galaxy. Combined with a determination of the neutral hydrogen column 
density via the Lya line, it is possible to determine the ISM abundances (e.g., [9-11]). We use S II, when available, 
as a measure of the metallicity, primarily since sulfur is not strongly depleted onto dust. We also place lower limits on 
the metallicity using Si II, Si IV and C II detections reported by Fynbo et al. [11]. 

Of the 18 GRBs in our sample, six have determined [S/H] values, five have no metallicity information, and six have 
upper limits on their metallicity from non-detections of S II as well as lower limits based on Si II or Si IV detections. For 
GRB 050908, the metallicity upper limit (Z < 10^ Zq) is not meaningful and we only report a lower limit based on a 
C II detection. Using these values we find a wide range of metallicities (using the Solar abundance values from Asplund 
et al. [21]): Z « 0.01 - 1.5 Zq for the Spitzer-dstected GRB hosts, which have stellar masses of 2 x 10^" Mq. We 
divide the GRBs with available metallicity information (either a metallicity detection or a bounded range) into two 
mass bins — the 3.6 /xm detections with ^2x10'" Mq (Group 1) and the objects included in the stack (Group 2). 
For each group, we perform a Monte Carlo simulation to estimate the mean metallicity, yielding (Zi) = —1.01 ±0.17 
and (Z2) = —1.52 ±0.12 (la intervals). For Group 1, the mean of the maximum inferred stellar masses is 3.7 x lO'*' 
Mq, while that of the masses inferred from the 70 Myr populations is 8.8 x 10^ Mq. To obtain mass estimates for 
Group 2, we scale our stack limit obtained for 1 1 non-detections by y/11/7. Using the mean maximum mass-to-light 
ratio of the objects in Group 2 yields an upper limit on the mean stellar mass of these seven objects of 3.7 x IO^Mq, 
while using the mean mass-to-light ratio at 70 Myr yields a mass limit of 9.4 x IO^Mq. 

In Figure 2, we present the absorption line metallicities plotted versus the stellar masses inferred from our Spitzer 
observations and the literature studies. The 6 GRB hosts that do not have any constraints on their ISM metallicity are 
excluded from this plot. We find that our two points fall below the observed relations at z < 3.5, providing evidence 
that the galaxy M*-Z relation continues to evolve at z ^ 3 — 5, with our stack range probing a somewhat lower mass 
scale than the LBG studies at z 3.1 - 3.5 [22, 23]. 




Stellar Mass (Solar units) 



FIGURE 2. Stellar mass plotted as a function of ISM metallicity for our sample and the 5 previously observed hosts 
(black=detections; gray=limits). ConfiiTned metallicities (and one metallicity lower limit) are indicated by filled symbols, while 
metallicity ranges are shown by dashed vertical lines with open symbols. The hatched regions designate 1 a and la intervals for 
estimates of the mean metallicity at two mass bins of ~ 2 x lO'" Mp (3.6 /xm detections) and < 3.7 x 10^ M.y (scaled stack limit - 
see text). The dashed vertical line indicates the upper limit on the mean mass of the stack for a 70 Myr population. These data are 
consistent with a decline in metallicity with lower stellar mass — an M*-Z relation. Also shown are the relations for z ~ 0.07 [3], 
z ~ 0.7 [4], z ~ 2.3 [5], and z 3. 1 — 3.5 [22, 23, filled squares]; the curves at z < 2.3 are M»-Z relations re-calibrated by Maiolino 
et al. [22]. Our two points at z 3 — 5 fall below these relations, suggesting that the M^-Z relation continues to evolve to z ~ 4. 



DISCUSSION AND CONCLUSIONS 



We present the first study of the galaxy mass-metallicity relation at redshifts of z ~ 3 — 5 using GRB afterglow 
absorption metalUcities and Spitzer follow-up observations. Five of the 20 GRB hosts in our sample are detected 
above a 3(J flux density threshold of 0.25 fiJy, corresponding to a typical stellar mass of ~ 2 X lO^o Mq. We further 
place a limit of < 3 x 10^ M© on the non-detected hosts based on a stacking analysis. 

Using the metalUcities of the host galaxies (inferred mainly from S ii) we find a range of ~ 10^ at a fixed stellar 
mass of-2x lO'^M.o.The mean metallicity at this mass scale is about 0.1 Zr. The mean metallicity associated with 
7 of the 1 1 the non-detected hosts, which have an upper limit of < 3.7 x lO'' M©, is Z < 0.03 Zq. Thus, there appears 
to be an overall decline in metallicity with decreasing mass, a hint of an M*-Z relation. Furthermore, our two points 
on the Mh,-Z relation lie below the relations at lower redshifts, suggesting that the relation continues to evolve at least 
to z 4. Clearly, additional observations are required to confirm and increase the statistical significance of this result. 
A sample of 20 additional GRBs at z > 3 from 2007 through the present is available for study. This will allow us to 
double the existing sample. 

The James Webb Space Telescope (JWST) will provide a much deeper view of the M*-Z relation at high redshift. 
For instance, the NIRCam instrument on JWST will allow us to detect GRB hosts at z ^ 3 down to a mass of 10** 
Mq, and at z ~ 6 to ~ 3 x 10* M© with similar integration times as our observations. Equally important, the NIRSpec 
instrument (1—5 ;Um) will allow us to determine emission-line metalUcities for some of the hosts, and hence to 
cross-calibrate the afterglow absorption metalUcities. With on-going afterglow absorption metalUcity measurements, 
the GRB sample will continue to play a key role in our study of high redshift galaxies. 
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